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We study both experimentally and theoretically modulation of light in a planar aligned deformed- 
helix ferroelectric liquid crystal (DHFLC) cell with subwavelength helix pitch, which is also known 
as a short-pitch DHFLC. In our experiments, azimuthal angle of the in-plane optical axis and 
electrically controlled parts of the principal in-plane refractive indices were measured as a function 
of voltage applied across the cell. Theoretical results giving the effective optical tensor of a short- 
pitch DHFLC expressed in terms of the smectic tilt angle and the refractive indices of FLC are used 
to ht the experimental data. Optical anisotropy of the FLC material is found to be weakly biaxial. 

For both the transmissive and rehective modes, the results of htting are applied to model phase 
and amplitude modulation of light in the DHFLC cell. We demonstrate that, if the thickness of the 
DHFLC layer is about 50 pm, the detrimental effect of held-induced rotation of the in-plane optical 
axes on the characteristics of an axicon designed using the DHFLC spatial light modulator in the 
reflective mode is negligible. 
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I. INTRODUCTION 

High-speed, low power consuming light modulation is 
in high demand for a variety of photonic devices used as 
building blocks of displays and optical information pro¬ 
cessors. These include tunable lenses, focusers, wavefront 
correctors and correlators M- 

For the most part of such devices, in addition to fast 
switching times, it is of crucial importance to have a 
27r modulation so that the phase can be smoothly tuned 
from zero to 27r. Liquid crystal (LC) spatial light mod¬ 
ulators (SLMs) are widely used as devices to modulate 
amplitude j)hase, or polarization of light waves in space 
and time Q- In LC-SLMs, nematic liquid crystals are 
among the most popular LC phases. However, nematics 
are known to have slow response time and, in addition, 
this slow response gets worse if the LC layer thickness 
increases in order to obtain the 27r phase modulation. 
Therefore, many efforts are in progress to optimize the 
various LC electro-optical modes for the high speed light 
modulation. 

Ferroelectric liquid crystals (FLCs) represent an alter¬ 
native and most promising chiral liquid crystal material 
which is characterized by very fast response time (a de¬ 
tailed descmtion of FLCs can be found, e.g., in mono¬ 
graphs [1, Q). Equilibrium orientational structures in 
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FLCs can be described as helical twisting patterns where 
FLC molecules align on average along a local unit direc¬ 
tor 

d = cos6>h-b sin0c, (1) 

where 9 is the smectic tilt angle; h is the twisting axis 
normal to the smectic layers and c T h is the c-director. 
The FLC director ([T|) lies on the smectic cone depicted 
in Fig. [T^ with the smectic tilt angle 9 and rotates in a 
helical fashion about a uniform twisting axis h forming 
the FLC helix with the helix pitch, P. This rotation is 
described by the azimuthal angle around the cone $ that 
specifies orientation of the c-director in the plane perpen¬ 
dicular to h and depends on the dimensionless coordinate 
along the twisting axis 

(j) = 2^(11-r)/P = qx, (2) 

where q = 2 tt/P is the helix twist wave number. 

Figure [1] illustrates the important case of a uniform 
lying FLC helix in the slab geometry with the smectic 
layers normal to the substrates and 

h = X, c = cos $ y-b sin $ z, 'E = Ez, (3) 

where E is the applied electric field which is linearly cou¬ 
pled to the spontaneous ferroelectric polarization 

Pg = PsP, p = h X c = cos $ z — sin $ y, (4) 

where p is the polarization unit vector. 

This is the geometry of surface stabilized FLCs (SS- 
FLCs) pioneered by Clark and Lagerwall in Ref. 0 
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Figure 1: (Color online) Geometry of (a) smectic cone 

where they studied electro-optic response of FLC cells 
confined between two parallel plates subject to homo¬ 
geneous boundary conditions and made thin enough to 
suppress the bulk FLC helix. Figure [Tb] also describes 
the geometry of deformed helix FLCs (DHFLCs) as it 
was originally introduced in [l^. The approach to light 
modulation that uses the electro-optical properties of he¬ 
lical structures in DHFLCs with subwavelength pitch also 
known as the short-pitch DHFLCs will be of our primary 
concern. 

In short-pitch DHFLC cells, the FLC helix is charac¬ 
terized by a submicron helix pitch, P < 1 /rm, and a 
relatively large tilt angle, 9 > 30°. Electro-optical re¬ 
sponse of DHFLC cells exhibits a number of peculiarities 
that make them useful for LC devices such as high speed 
spatial light modulators , colour-sequential liquid 

crystal display cells [la and optic fiber sensors [3 Il3 ■ 
The effects caused by electric-field-induced distortions of 
the helical structure underline the mode of operation of 
such cells. In a typical experimental setup, these effects 
are probed by performing measurements of the trans¬ 
mittance of normally incident light through a cell placed 
between crossed polarizers. 

In this article, our goal is to examine light modulation 
in planar aligned short-pitch DHFLC (PADHFLC) cells 
with uniform lying FLC helix (the twisting axis is par¬ 
allel to the substrates) and the related physical charac¬ 
teristics. This is the case that was studied theoretically 
in Refs [13, HH where the transfer matrix approach to 
polarization gratings was employed to define the effec¬ 
tive dielectric tensor of short-pitch DHFLCs. In par¬ 
ticular, it was found that, by contrast to the case of 
vertically aligned DHFLCs [T3, [i3|, the in-plane opti¬ 
cal axes of PADHFLCs sweep in the plane of the cell 
under the action of applied electric field thus producing 
changes in the polarization state of the incident light. 
More generally, biaxial anisotropy and rotation of the 
optical axes induced by the electric field in short-pitch 
DHFLC cells can be interpreted as the orientational Kerr 
effect [ilHllH. 

For detailed experimental characterization of this ef- 



and (b) planar aligned FLC cell with uniform lying helix. 

feet, we employ the experimental technique based on the 
Mach-Zehnder interferometer that goes beyond the lim¬ 
itations of the above mentioned standard experimental 
procedure and provides additional information on the 
principal refractive indices. Then we use the results 
as input parameters to study light modulation in the 
DHFLC cell operating in both the transmissive and re¬ 
flective modes. Our investigation into the effects of am¬ 
plitude modulation is based on the results of modelling of 
an axicon designed using the DHFLC spatial light mod¬ 
ulator in the reflective mode. 

An axicon as a cylindrically symmetric optical ele¬ 
ment that transforms an incident plane wave into a 
narrow beam of light along the optical axis have a 
long history dating back to the original papers by 
John N. McLeod [^, [23|. Laser beams propagating 
through axicons have two significant properties: (a) they 
generate a line focus, where the on-axis intensity stays 
high over much longer distances compared to focusing 
by conventional lenses; and (b) they generate ring-like 
intensity profiles in the far field. Both of these prop¬ 
erties proved useful in many applications such as atom 
guiding and trapping [^ . annular focusing in laser ma¬ 
chining 1^. g eneration of quasi-non-diffracting Bessel- 
like beams [264^ . subdiffraction limit imaging [ 23 | and 
optical micromanipulation [sol - l^ . 

The layout of the paper is as follows. In Sec. |lT] we 
introduce our notations, describe the recent theoretical 
results [2l| on the effective optical tensor of short-pitch 
biaxial FLCs, and discuss the orientational Kerr effect. 
Experimental details are given in Sec. Ell where we de¬ 
scribe the samples and the setup employed to perform 
measurements. The experimental data are fitted using 
the expression for the effective optical tensor. Sec. m 
deals with modulation of light in the DHFLC cells in 
both transmissive and reflective modes. Experimental re¬ 
sults are used to compute the amplitudes and phases of 
the components of transmitted and reflected light waves. 
DHFLC spatial light modulator acting as an axicon is 
modelled as a two-dimensional (2D) array of pixels which 
are DHFLC cells. In order to study the effects of am- 
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plitude modulation, the intensity distribution in the fo¬ 
cal plane of the axicon is evaluated for both DHFLC 
and ideal (no amplitude modulation) axicons. Finally, 
in Sec. we draw the results together and make some 
concluding remarks. Technical results are relegated to 
Appendices. 


II. ORIENTATIONAL KERR EFFECT 


In this section, we introduce notations and briefly dis¬ 
cuss the electro-optical properties of short-pitch DHFLC 
cells described by the effective dielectric (optical) tensor, 
£eff, defined in terms of averages over distorted FLC he¬ 
lical structure (see Appendix El- For the geometry 
of uniform lying FLC helix (see Fig. we recapitulate 
the analytical results for the optical tensor of a biaxial 
ferroelectric liquid crystal with subwavelength pitch . 
In the subsequent section, these results will be used to 
interpret the experimental data. 

We consider a FLC layer of thickness D with the z axis 
which, as is indicated in Fig.[TJ is normal to the bounding 
surfaces: z = 0 and z = D, and introduce the effective 
dielectric tensor, Seg, describing a homogenized DHFLC 
helical structure. 

For a biaxial FLC, the components of the dielectric 
tensor, e, are given by 

eij = e±6ij + (ei - ej.) didj + (e 2 - ej .)piPj 

= e±(<Jy + Uididj + U2PiPj), (5) 

where i,j S {a;,j/, z}, 5ij is the Kronecker symbol; di 
{pi) is the ith component of the FLC director (unit 
polarization vector) given by Eq. ((T|) (Eq. (|4])); m = 
(cj —e_L)/e_L = Aei/ej_ = Vi — l are the anisotropy parame¬ 
ters and ri = eije± (r 2 = e 2 /e_L) is the anisotropy (biaxi- 
ality) ratio. Note that, in the case of uniaxial anisotropy 
with it 2 = 0, the principal values of the dielectric ten¬ 
sor are: €2 = e_L and ei = e||, where n±_ = y/]de± 
(n|| = y^/Jel") is the ordinary (extraordinary) refractive 
index and the magnetic tensor of FLC is assumed to be 
isotropic with the magnetic permittivity p. We shall also 
assume that the medium surrounding the layer is opti¬ 
cally isotropic and is characterized by the dielectric con¬ 
stant Cm, the magnetic permittivity p^n and the refractive 
index Um = y/Jp^- 

Assuming that the pitch-to-wavelength ratio P/A is 
sufficiently small, the effective dielectric tensor can be 
expressed in terms of the averages over the pitch of dis¬ 
torted FLC helical structure. Explicit formulas for the 
components of the tensor are given in AppendixEl These 
formulas can be used to derive the effective optical tensor 
of homogenized short-pitch DHFLC cell for both verti¬ 
cally and planar aligned FLC helix [T^ [2lj . 

We concentrate on the geometry of planar aligned 
DHFLC helix shown in Fig. M For this geometry, the 
effective dielectric tensor can be written in the following 


form [2Ij : 

^xy^E 

^eff ~ I ^xy^E^ Cp -\- ^yyOi}^ 


(6) 

0 0 Cp — "lyyOi^ I 

The zero-field dielectric constants, Ch and Cp, that enter 
the tensor ([6]) are given by 

e?i/e_L = {nh/nxY = 


Ml cos^ 0 


r2 - 1 


,- 1/2 


+ y/r^ 

ep/e_L = {up/njfif = y/r^, 

M = r2 (m -b 1) = Ml sin^ 0-1-1. 


(7a) 

(7b) 
(7c) 

A similar result for the coupling coefficients jxx, lyy and 
7 a,p reads 


Zyjr 2 u a ■ a\i 

7xx/e_L = , r- , —;=^(mi cos0sm0) , 

(VM + y/r2Y 

7ru/f± = cos^sin^. 

Vu + y/ri^ 


(8a) 

(8b) 

(8c) 


Note that, following Ref. [l5|, we have introduced the 
electric field parameter 


o-E =XeE/Ps, (9) 

where xe = d{Pz)/dE is the dielectric susceptibility of 
the Goldstone mode [H, Is^ and P^ = Pg cos $. 

Note that the simplest averaging procedure previously 
used in Refs. [H, [T^ [l^ heavily relies on the first- 
order approximation where the director distortions are 
described by the term linearly proportional to the electric 
field. Quantitatively, the difficulty with this approach is 
that the linear approximation may not be suffice for accu¬ 
rate computing of the second-order contributions to the 
diagonal elements of the dielectric tensor ([S]) . In Ref. [2l| , 
the results (l6]l- (l8B were derived by using the modified 
averaging technique that allows high-order corrections to 
the dielectric tensor to be accurately estimated and im¬ 
proves agreement between the theory and the experimen¬ 
tal data in the high-held region. 

The dielectric tensor ([S|) is characterized by the three 
generally different principal values (eigenvalues) and the 
corresponding optical axes (eigenvectors) as follows 


Eeff = CzZ <8 Z-b e+d+0 d+-b e_d_ (g) d_, (10) 

£2 = = e^z ^ “ £p ~ lyy^E^ (H) 

= n| = e± -b [ixyCtEY ( 12 ) 

where 

£ = Y'^xP + £yy*^)/2 = £0 + [ixx + lyy)o?E/‘2., (13) 

A£ = - £^p®^)/2 = Aeo + ffixx - Xyy)a%l2, (14) 

£o = (£/t + £p)/2, Aeo = (£/t ~ £p)/2. (15) 
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The in-plane optical axes are given by 

d+ = cos'i/'dx-I-sin'0dy, d_ = z x d+, (16) 

2V'd = arg[Ae i'^xyOiE] ■ (17) 

From Eq. ([6]), it is clear that the zero-field dielectric 
tensor is uniaxially anisotropic with the optical axis di¬ 
rected along the twisting axis h = x. The applied electric 
field changes the principal values (see Eqs. (ITTl) and dT^ l 
so that the electric-field-induced anisotropy is generally 
biaxial. In addition, the in-plane principal optical axes 
are rotated about the vector of electric field, E || z, by 
the angle ipd given in Eq. (HU). 

In the low electric field region, the electrically induced 
part of the principal values is typically dominated by the 
Kerr-like nonlinear terms proportional to E'^ 

whereas the electric field dependence of the angle 4’d is 
approximately linear: E. This effect is caused by 

the electrically induced distortions of the helical structure 
and bears some resemblance to the electro-optic Kerr ef¬ 
fect. Following Refs. [Ulii , it will be referred to as the 
orientational Kerr ejfeet. 

It should be emphasized that this effect differs from 
the well-known Kerr effect which is a quadratic electro¬ 
optic effect related to the electrically induced birefrin¬ 
gence in optically isotropic (and transparent) materi¬ 
als [13. Similar to polymer stabilized blue phase liquid 
crystals [Mill, it is governed by the effective dielectric 
tensor of a nanostructured chiral smectic liquid crystal. 
This tensor is defined through averaging over the FLC 
orientational structure (see Appendix fSTl . 

III. EXPERIMENT 

In this section we present the experimental results on 
the principal refractive indices and orientation of the in¬ 
plane optical axis measured as a function of the applied 
electric field in DHFLC cells. 


A. Material and sample preparation 

In our experiments we used the FLC mixture FLC- 
587 (from P. N. Lebedev Physical Institute of Russian 
Academy of Sciences) as a material for the DHFLC 
layer. The FLC-587 is an eutectic mixture of the three 
compounds whose chemical structures are described in 
Ref. [ 13 . The phase transitions sequence of this FLC 
during heating up from the solid crystalline phase is: 

Cr+i^SmC* +H2!£,SmA* +i^Is, while cooling 


from smectic C* phase crystallization occurs around - 
10°C -15°C. The spontaneous polarization, Pg, and 

the helix pitch, P, at room temperature (22°C) are 
150 nC/cm^ and 150 nm, respectively. 

The cell was sandwiched between two glass substrates 
covered by ITO (indium tin oxide) and aligning films of 
the thickness 20 nm, and the gap was fixed by spacers 
at D = 50 mm. Geometry of the cells is schematically 
depicted in FigfTbl 

High quality planar alignment yielding the contrast ra¬ 
tio up to 200 : 1 was achieved using 4,4’-oxydianiline di¬ 
anhydride (PMDA-ODA) as aligning layers. Chemical 
formula of this polyimide after imidization is shown in 
Fig.m 



Figure 2: Chemical structure of the compound forming 
the aligning layer. 

PMDA-ODA dissolved in dimethyl-formamide (the 
concentration was about 0.2 wt%) was spin-coated onto 
the ITO surface. The polyimide film then was dried on 
the ITO substrate for 30-40 min at temperature 180°C, 
and subsequent imidization was done at temperature 
within the interval 275°C-290°C for about 1 h. 

Following the method of Ref. ^3 1 after cooling down, 
the polyimide films were rubbed with a cotton shred to 
provide aligning layers anisotropy. The FLC mixture was 
then injected into the cells in the isotropic phase by cap¬ 
illary action. 

Our task was to obtain regular helix alignment in the 
cell with the helix axis parallel to the glass plates. For 
this purpose, the FLC cell was subjected to an additional 
electrical training with a square-wave function of maxi¬ 
mum field amplitude ranged from 5 V/^m to 9 V//.tm and 
the frequency in the range between 0.5 Hz and 2 kHz [^ . 
The obtained alignment was inspected by observing tex¬ 
tures within the cell in a polarizing microscope. 


B. Measurement of azimuthal angle of in-plane 
optical axis 

In our experiments, we have used a low power He-Ne 
laser (A = 632.8 nm) as a light source. Initially, with¬ 
out applied voltage, the FLC cell was placed between 
the crossed polarizers and rotated so as to minimize the 
transmission of normally incident light. Then the cell was 
subjected to the time varying voltage of the symmetric 
square-wave form with the frequency of 40 Hz and the 
amplitude ranged from zero to 100 V. Under the action 
of the applied electric field, the in-plane optical axis ro- 
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Figure 3: (Color online) Azimuthal angle ipd of the major in-plane optical axis as a function of electric field for the 
DHFLC cell of thickness 13 = 50 /im filled with the FLC mixture FLC-587 [l^. The experimental points are marked 
by squares. Solid line represents the theoretical curve computed from Eq. (froi) with Ps/xe ~ 4.83 V//rm. 
Parameters of the mixture are: nj_ Ri 1.52 (e_L Ri 2.3) is the ordinary refractive index; Ri 1.77 (e|| = ei Ri 3.13) is 
the extraordinary refractive index, 9 = 35.5 deg is the tilt angle and r 2 = 1.05 is the biaxiality ratio. 


tates about the normal to the substrates (the z axis) and 
its azimuthal angle '0d changes. 

The angle ipd characterizing the electric-field-induced 
in-plane reorientation of the optical axis was measured 
by rotating the cell around the z axis and detecting the 
angle where the intensity of the transmitted light at pos¬ 
itive voltages is minimal. The experimental results for 
this angle obtained at different values of the voltage am¬ 
plitude are presented in Fig. [31 


C. Measurement of principal in-plane refractive 
indices 

In order to perform measurements of the principal val¬ 
ues of the in-plane refractive indices, n+ and n_ (see 
Eq. (HU)), we have used the well-known experimental 
method which is based on a Mach-Zehnder two-arm inter¬ 
ferometer (it is detailed in many textbooks such as [i^). 

In this method, as it can be seen in Fig. |4l a beam 
splitter divides a linear polarized incident light passed 
through the input polarizer into two paths and the FLC 
cell is placed in the path of the sample beam. The 
sample and reference beams are then recombined and 
pass through the output polarizer so that the interfering 
beams after the polarizer are collected by a photodiode. 

Given the amplitude of the voltage and the correspond¬ 
ing value of the principal axis azimuthal angle ipd-, the 
polarizers were rotated so that the polarization vector 
of the incident light is either parallel or perpendicular 
to the optical axis. In both cases, measurements giving 
the electrically controlled part of the corresponding re¬ 
fractive index were performed during the half-period of 


positive applied voltage. The results are shown in Fig. 0 


D. Results 

There are three optical characteristics of the DHFLC 
cell that we have measured in our experiments: the az¬ 
imuthal angle of the optical axis tpd and the electri¬ 
cally controlled parts of two principal refractive indices, 
An± = n±{E) — n±(0). The experimental data for elec¬ 
tric field dependence of the principal axis angle and An± 
are presented in FigjSjand Fig. [SJ respectively. 

We can now use formulas for n± [see Eq. m] and 
tpd [see Eq. (II3] to fit the experimental data. For this 
purpose, we assume that the FLC mixture is character¬ 
ized by the parameters: e_L r; 2.3 {n± = Uo ^ 1-52), 
e|| = 3.13 (n|| = Ue ~ 1.77) and 9 = 35.5°. The 
theoretical curves shown in Figs. |5]|ni are computed at 
PsIxe ~ 4.83 V/^m and r 2 = e 2 /e_L = 1-05. Interest¬ 
ingly, the value of the biaxiality ratio differs from unity 
and thus optical anisotropy of the mixture appears to be 
weakly biaxial. 

Note that the fitted values of the zero-field refractive 
indices are: Uh ~ 1.67 and Up r; 1.57. Figure [5] shows 
how the principal refractive indices, n± and n^, change 
with applied electric field. It is clear that electric-field- 
induced optical anisotropy is weakly biaxial. 


IV. LIGHT MODULATION IN DHFLC CELLS 

In this section, modulation of light in the DHFLC cells 
studied in the previous section will be of our primary con- 



Refractive index 
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Figure 4: Experimental setup for measurements of electrically controlled parts of the in-plane refractive indices of 

the DHFLC cell. 




(a) An_|_ = n_|_(iJ) — n_|_(0) (b) —An_ = ri—(0) — n—(E) 

Figure 5: (Color online) Electrically controlled differences of the principal in-plane refractive indices (a) An+ and 

(b) —An- as a function of applied electric field. 



1.57 


1.565 


X 

o 

■D 


.> 1.56 


a> 

OC 


1.555 


1.55, 



0.5 1 

Eiectric fieid, E (V/|rm) 

(b) 71— and rij^ 


Figure 6: (Color online) Electric field dependence of principal refractive indices. 


cern. For both the transmissive and reflective modes, the DHFLC modulator is found to be affected by the pres- 
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ence of amplitude modulation. We study how this mod¬ 
ulation influences the transformation characteristics of a 
DHFLC spatial light modulator operating as an axicon 
producing a ring-shaped far field distribution of light. 


A. Amplitudes and Phases 


Typically, in experiments dealing with the electro-optic 
response of DHFLC cells, the transmittance of normally 
incident light passing through crossed polarizers is mea¬ 
sured as a function of the applied electric field. In the 
case of normal incidence, the transmission and reflec¬ 
tion matrices can be easily obtained from the results of 
Refs. 111,1111,113 in the limit of the wave vectors with van¬ 
ishing tangential component, kp = 0. For our purposes, 
we shall need to write down the resultant expression for 
the transmission matrix 


_ n±jn- 

/m m 

where 




(19) 

( 20 ) 
( 21 ) 

(22a) 

(22b) 


and h = k^a_cD is the thickness parameter. Equation m 
defines the transmission matrix linking the vector ampli¬ 
tudes of incident and transmitted waves, Eq and Etransm, 
through the standard input-output relation 


Etransm — T(V'd)Eo, E. 


transm — 


,(transm) ' 
x 

,(transm) 

v y 


(23) 


When the incident wave is linearly polarized along the 
X axis (the helix axis) with 

Eo = , (24) 


the components of the transmitted wave can be written 
in the following form 

exp{i(j)^) = t+- - t_) sin^ ijji, 

Ei*~)/Eo = Ay expiicl^y) = sin(2V^d), (25) 

where A^^y and 4>x,y are the normalized amplitude and 
phase of the corresponding complex amplitude compo¬ 
nent, respectively. Then the transmittance coefficient de¬ 
scribing the intensity of the light passing through crossed 
polarizers is given by 






sin^(2V^d)- 


(26) 


Note that, under certain conditions such as |p±| ^ 1; 
both the transmission coefficients (EHl) and the transmit¬ 
tance (12611 can be approximated by simpler formulas 


t± Ri exp(f$±), = n±h, (27a) 

^ r; sin(A$/2) exp[z($+ -I- $_ -I- 7r)/2], (27b) 

\Ay\^ R! sin2(A$/2) sin2(2V>d), (27c) 


where A<I> = $+ — $_ is the difference in optical path 
of the ordinary and extraordinary waves known as the 
phase retardation. 

The complex valued components of the transmitted 
plane wave are thus characterized by the amplitudes and 
phases, A^^y and (t>x,y, given in Eq. (1^51) . These can now 
be computed as a function of the applied electric field 
by using our experimental data combined with the re¬ 
sults of fitting presented in Sec. IIII Dl Referring to Fig.[7l 
rotation of the optical axis combined with electric-field- 
induced change in the phase retardation produces varia¬ 
tions of the amplitudes Ax and Ay with electric field. 

For the component parallel to the helix axis. Ex, the 
amplitude modulation, however, has negligibly small ef¬ 
fect on the electric field dependence of the phase (j)x- 
It turned out that, despite amplitude modulation, this 
phase is close to <j)x ~ 4)-|.. 

From Eq. (l27cD . it might be concluded that the phase 
of the component perpendicular to the helix axis, (jiy is 
given by ($+ -|- <!)_ -|- 7r)/2. It should be noted that, 
in the approximation described by Eq. (j27al) . when the 
factor sin(A$/2) changes its sign at the points where 
the amplitude Ay (and Txy) is zero, the phase should 
experience jump-like behaviour with A^j, = ±7r. These 
jumps are not shown in Fig. [7] as, in real experiments, 
the amplitude Ay never reaches zero due to scattering 
effects. 

We have also computed the amplitudes and phases of 
the reflected wave for the DHFLC cell operating in the 
reflective mode with the mirror placed at the exit face 
of the cell. It can be shown that, for the case of normal 
incidence, the reflection matrix can be written as follows 
(details are relegated to Appendix IbI) 

R-(V’ci) = ^ era -h ^ Rt(-2V’d) (28) 

^ P±+P± exp(2f$±) ^ Rr-P± , , 

1 + P±P± exp(2f$±) ’ 1 — RrP± ’ 


where Rr is the reflection coefficient of the mirror given 

by 




1 — exp{2inrhr) 

1 — Px exp(2inrhr) ’ 


Tlj./ Pr ^m/ pm 
rir /pr n^a/Pm 

(30) 


Ur (pr) and hr = kvacDr (Dr) is the mirror refractive in¬ 
dex (magnetic permittivity) and the thickness parameter 
(thickness) of the mirror, respectively. 
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Figure 7: (Color online) Amplitudes, and Ay and phases, (px and (j>y, for the components of the transmitted light 
vector amplitude Etransm/Ao computed as a function of applied electric field in the transmissive mode. 




Figure 8: (Color online) Amplitudes, Ax and Ay and phases, px and py, for the components of the reflected light 
vector amplitude Ej-efl/Eo computed as a function of applied electric field in the reflective mode. Dotted lines 
represent the curves computed from Eq. (EH) with the doubled layer thickness {D is replaced by 2D). Solid curves 

are evaluated using the reflection matrix (E51) and Eq. (15^ . 


In the reflective mode, for the incident wave (12411 . the input-output relation 


(refl) 

(refl) 


Erefl — I^(V^d)Eo, Efefl — 


(31) 
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gives the components of the reflected wave 


alent to the phase profile of an axicon given by 


exp(i(/)^) = r+ - (r+ - r_) sin^ ipd, 

= Ay exp{i(j)y) = ^ sin(2V)d)- (32) 

The curves depicted in Fig. |S] as solid lines are com¬ 
puted for the silver mirror of the thickness Dr = 0.15 /im 
which is characterized by the complex valued refractive 
index [i^ Ur ~ 0.16-1-3.at the wavelength A = 633 nm. 
The magnitude and the phase of the reflection coeffi¬ 
cient O, Rr = \Rr\exp{i(f>r), then can be estimated 
at about |i?r| ~ 0.97 and 4>r ~ 43°. 

Alternatively, the reflective mode can be described in 
the double-layer approximation using the transmission 
matrix (HU) where the thickness of the cell is doubled 
and D is replaced by 2D. The results computed in this 
approximation are shown in Fig.|5]as dotted lines. Inter¬ 
estingly, for the amplitude modulation of the component 
Ay, the exact and approximated results are in excellent 
agreement. By contrast, the curves for the magnitude of 
the component A^ [see Fig. |Sl[a)] noticeably differ from 
each other. It comes as no surprise that the double-layer 
approximation, which is based on the assumption of per¬ 
fect reflection, overestimates \Ax\. 

Similar to the transmissive mode, the results for the 
phases, (px and (f>y, can be easily understood in the lim¬ 
iting case where |p±| are small and the reflection coeffi¬ 
cients r± defined in Eq. (1^ are close to exp(2i$±). 
From Eq. (15^ . it can be inferred that (j)x ~ 2^+ + (pr 
and (j)y ~ $+-|-4>_-|-^!)r- + 37r/2, whereas the double¬ 
layer approximation gives the relations: (px ~ 2$-|_ and 
py Ri $+-|-4>_-|-7r/2. The curves plotted in Figs.[51[c)-(d) 
clearly show the phase shift introduced by the mirror and 
the minus sign on the right hand side of equation for 
the component Ay. It can also be seen that phase modu¬ 
lation is about two times larger as compared to the case 
of the transmissive mode. So, the electric field required 
to reach 2tt modulation is reduced by the factor about 
1.4. 


B. Axicon 

The liquid crystal spatial light modulators are exten¬ 
sively used for formation of light wavefields with a speci¬ 
fied spatial distribution of intensity. The above discussed 
effects of amplitude modulation may affect quality of the 
resulting light field. In order to estimate how these ef¬ 
fects influence a reflective DHFLC modulator, we have 
modelled the DHELC-SLM as a 2D array of 256 x 256 
pixels each of the area 200 x 200 /rm^. It should be noted 
that in our modelling the electric field is assumed to be 
uniform. In our case, the pixel size to the cell thickness 
ratio is large enough for this assumption to be justified. 

In our calculations, we assumed that the incident wave 
is linearly polarized along the helix axis (x axis) and the 
phase distribution in the plane of the modulator is equiv- 


— 4?axicon 


PO 

i?’ 


Po 



(33) 


where k = 2 tt/\ is the wavenumber; / is the focal length; 
R is the axicon radius and 4>axicon is the phase modula¬ 
tion depth of the axicon. The x component of the light 
field in the focal plane of lens can be computed from the 
Fresnel-Kirchhoff diffraction formula taken in the far field 
(Fraunhofer) approximation (see, e.g.. Ref. 0 ) 


Ex{x,y) = 
X exp 


ex 


:p[ifc/+ i-oo 


iXf 

ik 


/ 7 

J —oo J —oo 


Ax{xo,yo) 


ipx(xo,yo) - -J {xxo + yyo} dxodj/o- (34) 


The amplitude and phase distributions, Ax{xo, yo) and 
Pxixo,yo) that enter the integrand on the right hand side 
of Eq. (IMl) are approximated as follows: (a) we recast 
the phase profile of an axicon characterized by the speci¬ 
fied depth of modulation into the step-like form with the 
step height equal to 27r; (b) the spatial phase distribu¬ 
tion is then discretized along the coordinates x and y 
with /S.X = Ay = 200 /rm so that each pixel is charac¬ 
terized by the constant phase equal to the value of px in 
its center; (c) for each value of the phase px, we com¬ 
pute the corresponding value of the amplitude Ax and 
derive the discretized distribution of the amplitude. The 
final step involves using the standard FFT (fast Fourier 
transform) technique [45| | to evaluate of the integral on 
the right hand side of Eq. (EH). 

Figure E] shows the phase and amplitude distributions 
that were computed for the modulation depth, 4>axicon, 
equal to 167r. The intensity distribution in the focal plane 
and the corresponding x-dependence of the intensity are 
presented in Fig. [TU] and in Fig. [TTl respectively. In 
Fig. El the case of DHFLC axicon (dashed line) (see 
Fig. llOall is compared with the curve (solid line) com¬ 
puted for an ideal axicon without amplitude modulation 
(see Fig. llObp . The latter implies that the amplitude 
-^x{xo,yo) is assumed to be constant. For the curves 
depicted in Fig. El loss of power caused by amplitude 
modulation can be estimated to be below 5%. 

Note that the results shown in Figs [5] - El are com¬ 
puted at the DHFLC cell thickness taken to be about 
50 /im. We have also found that, when the thickness of 
the DHFLC cells is halved, reducing the thickness results 
in an increase of the power loss that can be estimated at 
about 11%. 


V. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied optical properties of 
the planar aligned DHFLC cells (geometry is shown in 
Fig. P) that govern modulation of light in such cells. 
In our experiments, we have combined the experimen¬ 
tal data for electric field dependence of in-plane optical 
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(a) Distribution of (px (b) Distribution of Ax 

Figure 9: Profiles of (a) phase, (fix, and (b) amplitude, Ax, in the transverse plane of the DHFLC spatial light 
modulator. Grayscale with 256 levels of gray encodes the values of (fx and Ax ranged from zero (black color) to 27r 

and unity (white color), respectively. 



(a) Intensity in focal plane of the DHFLC axicon 


(b) Intensity in focal plane of the ideal axicon 


Figure 10: Grayscale representation of the intensity distribution in the focal plane of (a) the DHFLG and (b) ideal 
(no amplitude modulation) axicons operating in the reflective mode. The cell thickness is about 50 /im. 


axis orientation characterized by the azimuthal angle ifd 
(see Fig. [3]) with the experimental technique based on 
the Mach-Zehndler interferometer to measure the prin¬ 
cipal refractive indices as a function of the electric field 
(see Fig. O. 

General expression for the effective optical tensor of a 
DHFLG cell derived in Ref. [2l| [see Eq. ([5])] was used to 
fit the experimental curves. The obtained parameters of 
the FLG mixture FLC-587 show that this material is a 


weakly biaxial FLG. Similarly, the resulting electric field 
dependence of the principal effective refractive indices of 
the DHFLC cell plotted in Fig. [5] clearly indicate electric- 
field-induced optical biaxiality of the effective dielectric 
tensor ®. Note that optical biaxiality of the smectic C* 
phase was previously reported in Refs 

We have analyzed modulation of light in the DHFLC 
cell under consideration by computing the amplitudes 
and phases of the components of light wave transmit- 
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Figure 11: (Color online) Intensity in the focal plane as a function of x for the DHFLC and ideal operating in the 
reflective mode. Solid line represents the distribution computed for the case of an ideal axicon without amplitude 

modulation. The cell thickness is about 50 ^m. 


ted through (reflected from) the cell in the transmissive 
(reflective) mode described using the transfer matrix ap¬ 
proach [^. It is found that, in addition to phase mod¬ 
ulation, electric-field-induced birefringence and rotation 
of the in-plane optical axes generally result in the effects 
of amplitude modulation. Additionally, we have calcu¬ 
lated the far field intensity distribution formed by the 
DHFLC spatial light modulator operating as an axicon. 
From a comparison between the results obtained for the 
ideal (no amplitude modulation) axicon and the DHFLC 
modulator, it can be concluded that the characteristics 
of the DHFLC axicon are very close to the ones of the 
ideal axicon. In particular, for the cell thickness about 
50 /im the power loss due to amplitude modulation is es¬ 
timated at about 5%, whereas it increases up to 11% at 
the halved thickness. For many applications, such values 
of the power loss can be regarded as acceptable. Note 
that, according to Ref. [i^, the vortex lightwave fields 
shaped in the form of certain curves, which are of impor¬ 
tance for efficient laser micromanipulation, can be gener¬ 
ated using the 2D array made up of, at least, 128 x 128 
amplitude-phase diffractive elements. 
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Appendix A: Effective dielectric tensor 


According to Ref. [^, the effective dielectric tensor 


^efF — 


/g(eff) 

^XX 

Jeff) 

tyx 


V 


Jeff) 

txy 

,(eff) 


(eff) (eff) 




Jeff)N 

^XZ 

,(eff) 

tyz 

,(eff) 

^ZZ 


(Al) 


can be expressed in terms of the averages 

Vzz = (ej^) = ej^([l -b uidl + U2pl]~^), 

P _ /, /, \ _ / Uldzda + U2PzPa\ 

Pza — \i-zal^zz) “ \ i , ,2 , 2 / 

\ 1 -b Uid% + '^2Pz / 


(A2) 

(A3) 


In addition, the typical response time of the DHFLC 
modulator based on FLC-587 is known to be around 
100 ps at the phase shift equal to tt da. So, it might be 
concluded that all essential prerequisites are in place for 
elaboration of the DHFLC spatial light modulator oper¬ 
ating almost as the ideal axicon at the phase modulation 
frequency about 1 kHz. 


where (...) = (.. .)^ = {2-k) ^ ...dcj) and a S {x, y}, 

as follows 

= Pza/Vzz, 

^ = (4^^) + I3zal3z0/Vzz, (A4) 
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where are the components of the averaged tensor, 

(£p): 


( 4 ?) = 


— eo ( <^a/3 + 


Uidgdp + U2PaPl3 + UiU2qaqp 
I-\-uid1 + U2 pI 
Qa = Pzda - d:,Pg, a,l3€{x,y}, 


(AS) 

(A6) 


describing effective in-plane anisotropy that governs 
propagation of normally incident plane waves. 

General formulas (IMI-dMI) give the zero-order ap¬ 
proximation for homogeneous models describing the op¬ 
tical properties of short-pitch DHFLCs [H, ■ 


Appendix B: Derivation of the reflection matrix 


In this Appendix, our task is to derive the reflection 
matrix (I28|l for the plane wave normally incident on the 
DHFLC cell which exit face is covered with a thin re¬ 
flecting layer (mirror). For this purpose, we shall use 
the transfer matrix approach in the form presented in 
Refs. l2l|,[53,[iil. We begin with the 4x4 transfer ma¬ 
trix 


W‘0=WPLOW„ (Bl) 


corresponding result for the reflection matrix reads 

RFLC=W^r"^^[Wil 

= (T3Rt(V’d) 

T] 

(FLC) 1 - exp{2in±h) 

r± - P±- 


( (FLC) 

.^(FLC) ) R-t( —-i/'d), 


0 


1 — exp(2m±^) 

For the reflecting layer, the results 
— rwG)i-i — ' 

1 - Pr 


(B6) 

(B7) 


T. = [WiV]-^ = 


T, = 


1 — Pr exp{2inrhr) 


exp{inrhr), 


(B8a) 

(B8b) 


R,, = w^v[w 


rG)i-i 




(B8c) 


Rr — Pr 


1 — exp(2inrhr) 

1 — p'^ exp{2inrhr)’ Ur/pr+nm/p 


r/Pr Pm/Pn 


(B8d) 


can be obtained from the above formulas for the FLC 
layer by replacing {n+, n_, h = fcvacTI} with {ur, hr = 
hvac Rr} ■ 

For a non-absorbing FLC material, the unitarity rela¬ 
tions imply the symmetry conditions 211. 


[W, 


(flc)it ^ 


[W 


(FLC)iT 
21 J 


= -W 


[W(f «)]t = <T3W(f 


[W, 


(FLC)nt _ 


= -^3W(f "^^^3, 


(B9) 

(BIO) 

(Bll) 


where W^- are 2x2 block-matrices, expressed as a prod¬ 
uct of the transfer matrices of the FLC cell and the mir¬ 
ror, WpLC and Wr, respectively. Given the transfer ma¬ 
trix (EH), the transmission and reflection matrices, T and 
R, can generally be computed from the formulas (see, 
e.g.. Ref. dH) 

T = wri\ R = W2 iW3A' (B2) 

that relate T and R with the block matrices Wn and 

W 21 . 

Applying Eq. (IB2I) to the case of light normally imping¬ 
ing upon the FLC layer characterized by the dielectric 
tensor m gives the transmission matrix 


where a dagger and the superscript T will denote Hermi- 
tian conjugation and matrix transposition, respectively. 
The conditions (IB9I) can now be used to express the 
block-matrix Wn in terms of the transmission and re¬ 
flection matrices as follows 

Wn = 

= -b 

= TFlc[l2-RFLcRr]T-l. (B12) 

Equation (IB12I1 combined with the relation (ESI) imme¬ 
diately give the transmission matrix in the following gen¬ 
eral form 


Tflc = 

C(FLC) Q \ 

= Rt('0d) ( ^(FLC) j Rf(~V'd), 


RFLC) _ 
— 


1 — p^ exp(2m±h) 


exp{in±h), 


P± 


n±lp- TlpYl /P m 
n±lp + 

Tlpn /p m 


T = T,[l2-Ri:LcRr]"^TFLC. (B13) 
Similarly, using the symmetry relations (IBIOII and (|Bllll . 


(B3) 

we can obtain the following result for the block-matrix 
W 21 : 

(B4) 

W 21 = wW -f w« 


(B5) 

= -b 

= 0'3[Tflc] ~ ['’’aRFLcJ^lTr 

(B14) 


where Rt('0d) is the rotation matrix [see Eq. (I22bll ]. Substituting formulas (IB12I) and (IB14I) into Eq. (IB2I1 

which is identical to the matrix given in Eq. (HU). The gives the resulting expression for the reflection matrix 
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in the following form: 

R = cr3[Tj;,Lc]”H'^3Rr - [o-sRfLc]'''} 

X [I 2 — Rpn^R-r] TpLC- (B15) 

Formulas (IBlSp and (|B15I) are quite general. In par¬ 
ticular, they can be applied to any non-absorbing nni- 
formly anisotropic layer represented by the FLC mate¬ 
rial. For the case of normal incidence, from Eqs. (1531) - 
(IB8|) . it follows that, similar to the matrices TFLc(V'd) 


and RFLc(V’d), the matrices T('0d) and R(V'd) can be 
recast into the following factorized form: 

T(V^d) = Rt(V’d)T(0)Rt(-V'd), (B16a) 

R(V'd) = Rt(-?/'d)R(0)Rt(-V’d)- (B16b) 

The matrices, T(0) and R(0), are diagonal and can be 
easily computed. After algebraic manipulations, it is not 
difficult to derive the analytical expression for the trans¬ 
mission matrix given in Eqs. (EH])-®- 
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